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Abstract 

We prove that the class of GUM matrices is the largest class of bi-potential 
matrices stable under Hadamard increasing functions. We also show that any power 
a > 1, in the sense of Hadamard functions, of an inverse M-matrix is also inverse M- 
matrix showing a conjecture stated in Neumann jEJ. We study the class of filtered 
matrices, which include naturally the GUM matrices, and present some sufficient 
conditions for a filtered matrix to be a bi-potential. 

AMS subject classification: 15A48, 15A51, 60J45. 
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1 Introduction and Basic Notations 

A nonnegative matrix U is said to be a potential if it is nonsingular and its inverse 
satisfies 

Ur 1 < for j, U^>0 

3 

that is U" 1 is an M-matrix which is row diagonally dominant. We denote this class 
of matrices by V. In addition we say that U is a bi-potential if U" 1 is also column 
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diagonally dominant. We denote this class by biV. We note that V, biV are contained in 
.M _1 the class of inverses of M-matrices. 

The class of potentials matrices play an important role in probability theory. They 
represent the potential (from where we have taken the name) of a transient continuous 
time Markov Chain (X t ) t >o, with generator — U~ l . That is, 

/»oo /*oo 

U iS = I (e- u ~ H ) ij dt= W i {X t = j}dt 
Jo Jo 

is the mean expected time expended at site j when starting the chain at site i. Clearly 
U is a bi-potential if both U and U' are potentials. 

To get a discrete time interpretation take K = maxjf/j" 1 }- For any k > K the matrix 

i 

Pfc = I — \U~ l is nonnegative, sub-stochastic and 

U- 1 =k(l-P k ), 

If we can take k = 1, then U~ l = I — P (with P — Pi) and U is the mean expected 
number of visits of a Markov chain (Y n ) n£ ^ whose transition probability is given by P, in 
fact 

n>0 n>0 

We notice that if U is a potential then it is diagonally dominant on each row, in the 
sense that for all i,j we have Uu > Uji. The probabilistic prove of this fact is that 

where fji < 1 is the probability that the Markov process {X t ), starting from j ever reaches 
the state i. If U is a bi-potential, then it is also diagonally dominant on each column. In 
this case we just say it is diagonally dominant. 

For any nonnegative matrix U we define the following quantity 

t(U) = inf > : I + tU £ biV}, 

which is invariant under permutations that is t(U) = t(UUU'). We point out that if U 
is a positive matrix then r(U) > 0. We shall study some properties of this function r. In 
particular we are interested in matrices for which r(U) = oo. The next result shows that 
t(U) = oo is a generalization of the class biV. 

Proposition 1.1 Assume U is a nonnegative matrix, which is nonsingular and r{U) = 
oo, then U G biV . 

Proof. It is direct from the observation that 

tii + tuy 1 -> U' 1 . 

t— >oo 

□ 
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Remark 1.1 We shall prove later on that the reciprocal is also true: ifU is in the class 
biP , then t(U) = oo. 

The following notion will play an important role in this article. 

Definition 1.1 Given a matrix B we say that a vector fi is a right equilibrium potential 



where 1 is the constant vector of ones. Similarly it is defined the notion of a left equilibrium 
potential, which is the right equilibrium potential for B' . When B is nonsingular we denote 
the unique right and left equilibrium potentials by /i# and v^. 

We denote by Ji = l'/i the total mass of \x. In the nonsingular case, it is not difficult 
to see that v = fx. 

Notice that for a matrix U G biV the right and left equilibrium potentials are nonneg- 
ative. This is exactly the same as the fact that the inverse is row and column diagonally 
dominant. 

Definition 1.2 Constant Block Form (CBF) matrices can be defined recursively in the 
following way: given two CBF matrices A, B of sizes p andn—p respectively, and numbers 
a,/3 we produce the new CBF matrix by 



where the vector l p is the vector of ones of size p. We also say that U is in increasing 
CBF if min{A,B} > min{a, (3} . 

We recall the following two definitions introduced in ^JJ and [Til, that generalize the 
concept of ultrametric matrices introduced in PU] (see also [T3]). 

Definition 1.3 A nonnegative CBF matrix U is in Nested Block Form (N~BF) if in 
A, B are NBF matrices and 

• < a < (3; 

• min^y, Aji} > a and mm{B M , B lk } > a; 

• m&xiAijjAji} > (3 and m&x{B kh B lk } > (3. 

Definition 1.4 A nonnegative matrix U of size n, is said to be Generalized Ultrametric 
Matrix (GUM) if it is diagonally dominant that is for all i,j it holds Ua > max{[Zy, Uji}, 
and n < 2, or n > 2 and every three distinct elements i,j,k has a preferred element. 
Assume this element is i, which means 



B/2=l 




(1.1) 
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• Uji Ufa, 

• min{Uj k ,U kj } > mm{Uji, Uij}; 

• max{U jk , U kj } > max{Uji, U^}. 

By definition the transpose of a GUM matrix is also GUM. We note that an ultrametric 
matrix is a symmetric GUM. The study of the chain associated to an ultrametric matrix 
was done in [5] and for a GUM in jH]. 

In the next Theorem we summarize the main results in and concerning GUM 
matrices. 

Theorem 1.1 Let U be a nonnegative matrix. 

• U is a GUM matrix iff it is permutation similar to a NBF. 

• If U is a GUM, then it is nonsingular iff it does not contain a row of zeros and no 
two rows are the same. 

• If U is a non singular GUM then U G biV. 

It is clear that if U is a GUM then I + tU is a nonsingular GUM. In particular r(U) = oo. 
We introduce a main object of this article. 

Definition 1.5 Given a function f and a matrix U, the matrix f(U) is defined as 
f(U)ij = f(Uij). We shall say that f(U) is a Hadamard function ofU. 

Given two matrices A, B of the same size we denote by A © B the Hadamard product 
of them. So (A B)^ = AijBij. 

Given a vector a we denote by D a the diagonal matrix, whose diagonal is a. For 
example we have D a D b = D a Q D b = D aQb . 

The class of CBF matrices (and its permutations) is closed under Hadamard functions. 
Similarly, the class of increasing CBF (and its permutations) is closed under increasing 
Hadamard functions. 

On the other hand the class of NBF, and therefore also the class of GUM matrices, 
is stable under Hadamard nonnegative increasing functions. We summarize this result in 
the following proposition. 

Proposition 1.2 Assume U is a GUM and f : IR + — > M + is an increasing function. Then 
f(U) is a GUM. In particular r(f(U)) = oo, and if f(U) is nonsingular then f (U) G biV . 
A sufficient condition for f{U) to be nonsingular is that U is nonsingular and f is strictly 
increasing. 
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Proof. It is clear that f(U) is a GUM matrix and therefore r(f(U)) = oo. Then, from 
Proposition II . II we have that f(U) £ biP as long as it is nonsingular. If U is nonsingular 
then it does not contain a row (or column) of zeros and there are not two equal rows 
(or columns). This condition is stable under strictly non-negative functions, so the result 
follows. □ 

One of our main results is a sort of reciprocal of the previous one. We shall prove that 
if r(f(U)) = oo for all increasing nonnegative functions /, then U must be a GUM (see 
Theorem in)). 

The last concept we need for our work is the following. 
Definition 1.6 We say that a nonnegative matrix U is in class T if 



and I + t(U) U is nonsingular whenever r(U) < oo. 

We shall prove that every nonnegative matrix U which is a permutation of an increasing 
CBF, is in class T . 

We remark here that our purpose is to study Hadamard functions of matrices and not 
spectral functions of matrices, which are quite different concepts. For spectral functions of 
matrices there are deep and beautiful results for the same classes of matrices we consider 
here. See for example the work of Bouleau [3] for filtered operators. For M matrices see 
the works of Varga [T^], Micchelli and Willoughby [T2], Ando pQ, Fiedler and Schneider 
jHj, and the recent work of Bapat, Catral and Neummann j2] for M-matrices and inverse 
M-matrices. 



Theorem 2.1 Assume U £ V and f : R + — > IR + is a nonnegative strictly increasing 
convex function. Then f(U) is nonsingular and det(f(U)) > 0. Also f(U) has a right 
nonnegative equilibrium potential. Moreover if /(0) = we have M = U~ 1 f(U) is an 
M-matrix. If U £ biV then f(U) also has a left nonnegative equilibrium potential. 

Note that H = /([/) _1 is not necessarily a Z-matrix, that is for some % ^ j it can 
happen that Hij > 0, as the following example will show. Therefore the existence of a 
nonnegative right equilibrium potential, which is 



t(U) 




2 Main Results 



VZ //„ • ^//; ; M). 
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does not imply that the inverse is row diagonally dominant, that is 

Vi Ha > ^2 \Hij\- 

Example 2.1 Consider the matrix 




Then U = (I — P) 1 E biV. Consider the strictly convex function f(x) = x 2 — cos(x) + 1. 
A numerical computation gives 

I 0.3590 -0.0975 0.0027\ 
(f(U)y 1 « -0.0975 0.2372 -0.0975 , 
\ 0.0027 -0.0975 0.3590/ 

which is not a Z-matrix. 

We denote by the Hadamard transformation of U under f(x) = x a . In particular 
=UQU. One of our main results is the stability of M. 1 under powers. This solves 
a conjecture stated on [T5] . 

Theorem 2.2 Assume U E M" 1 and a > 1. Then E M~ l . If U" 1 E V then 
([/(«))-! eV. If U E biV then E biV. 

The previous result has the following probabilistic interpretation. If U is the potential 
of a transient continuous time Markov process then is also the potential of a transient 
continuous time Markov process. In the next result we show the same is true for a 
potential of a Markov chain. An interesting open question is what is the relation between 
the Markov chain associated to U and the one associated to U^ a \ 

Theorem 2.3 Assume that U^ 1 = I — P where P is a submarkov kernel, that is P > 0, 
PI < 1. Then for all a > 1 there is a submarkov kernel Q(a) such that (U^)^ 1 = 
I - Q(a). Moreover if P'l < 1 then Q(a)'l < 1. 

The next result establishes that the class of GUM matrices is the largest class of 
potentials stable under increasing Hadamard functions. 

Theorem 2.4 Let U be a nonnegative matrix such that r(f(U)) = oo, for all increasing 
nonnegative functions f. Then, U must be a GUM. 
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Example 2.2 Given a, b,c,d E M + consider the non-singular matrix 



(1 








o\ 





1 








a 


b 


1 







d 





V 



For all increasing non-negative functions f and all t > 0: (I + £/({7)) -1 is an M -matrix, 
while U is not a GUM. Moreover, U is not a permutation of an increasing CBF. This 
shows that the last Theorem does not hold if we replace the class biV by the class M~ l . 

Theorem 2.5 Let U G biV and f : IR + — > IR + be a strictly increasing convex function. 
f(U) is in biV if and only if f(U) belongs to class T. 

Theorem 2.6 If U is a nonnegative increasing CBF matrix then U is in class T. 



As a corollary of the two previous theorems we obtain the following important result. 

Theorem 2.7 Assume that U G biV is an increasing CBF matrix, and f : IR + — > R + is 
a nonnegative strictly increasing convex function. Then f(U) G biV. 



3 Proof of Theorems ED, E21 EH and ET5 

Let us start with a useful lemma. 

Lemma 3.1 Assume U G Mr 1 . Then for all t > we have (I + tU) G M~ l . Also if 
U G V so is (I + tU) and its right equilibrium potential is strictly positive. In particular 
if U G biV then so is I + tU and its equilibrium potentials are strictly positive. Similarly, 
let < s < t and assume I + tUE biV, then I + sll G biV and its equilibrium potentials 
are strictly positive. 

Proof. For some k > large enough, U^ 1 = k(I — N) where N >0 (and Nl < 1 for the 
row diagonally dominant case). In what follows we can assume that k — 1 (it is enough 
to consider the matrix kll instead of U). 

From the equality (I - N)(l + iV + TV 2 H N p ) — I — N p+1 we get that 

I + N + N 2 + ■ ■ ■ N p = U(I — N p+1 ) < U, 

oo 

and we deduce that the series ^2 N l is convergent and its limit is U. 

i=i 

Consider now the matrix 

w = (III __l_ JV )- 1 _ I[ ) =i g(_^Viv' 



We have that N t > (and N t 1 < 1 whenever Nl < 1). Therefore the matrix I — Nt is an 
M- matrix (which is row diagonally dominant when M is so) . On the other hand we have 

I + tU = I + t(I - iV)" 1 = (tl + 1 — AT) (I — N)' 1 = (1 + f) MI - j^n) (I - AO" 1 , 

from where we deduce that I + tU is nonsingular and its inverse is 
(I + (£/)-■ =T i.(I_J V )(l_ T i.JV)-' 



^((I-T^r'-^I-T^r 1 ) 

(oo oo \ 

j2( 1 + t) -i N i_j2 {1 + t) -i N i + i\ 
1=0 1=0 J 



The only thing left to prove is that A^l < 1 in the row diagonally dominant case, that is 
when Nl < 1. For that it is enough to prove that iV'l < 1 for large /. From the equality 

oo oo 

U = Y^N l , we deduce that ^2 N l l < oo, and therefore iV'l tend to zero as / — ► oo. This 

i=i i=i 
proves the claim. 

oo 

When k is not 1 we have the following equality (I + tffp 1 = ^(1 - f E(^) ZiV '), 

i=i 

where iV = I — \U~ l . 

Finally, assume that I + tU E biV. Hence l + /3(I+ tU) E biV for all j3 > 0. This 
implies that 

I + -^—t U E biV. 
1 + P 

Now, it is enough to take (3 > such that s = j+gt. □ 

This Lemma has two immediate important consequences. 

Corollary 3.1 If U E biV then t(U) = oo. 

Corollary 3.2 Let U a nonnegative matrix, then 

t(U) = sup{t > : I + tU E biV} 

Proof. It is clear that r{U) < sup{t > : I + tU E biV}. On the other hand if 
I + tU E biV then we get I + sll E biV for all < s < t. This fact and the definition of 
t(U), implies the result. □ 

Proof. (Theorem EHJ) We first assume that /(0) = 0. We have that U' 1 = k(I - P), 

for some k > and P a sub-stochastic matrix. Without loss of generality we can assume 
k = 1, because it is enough to consider kU instead of U and f(x) = f(x/k) instead of /. 



S 



Consider M = (U 1 f(U)). Take i ^ j and compute 

Ma = iU~ l f{U)) l3 = (1 -PiJfiUij) - 5> ifc /(E/iy)- 

Since 1 — pa — ^ Pifc > 0, which is equivalent to ^Pik < 1, and / is convex we obtain 

kj^i k 



5>* ) /(°) + > / ( ) 

k J k \ k / 



f{Ua). 



The last equality follows from the fact that U — I — P. This shows that < 0. 
Consider now a positive vector a; such that y' = x'U' 1 > (see [H], Theorem 2.5.3). Then 

x'M = x'U^fiU) = y'f(U) > 0, 

which implies, by the same cited theorem, that M is a M-matrix. In particular M is 
nonsingular and det(M) > 0. So, f(U) is nonsingular and det(f(U)) > 0. Consider now 
p the right equilibrium potential of f(U). We have 

Mp = U~ l f(U)p = U~H = pu > 0, 

then p = M~ x pu > 0, because M _1 is a nonnegative matrix. This means that f(U) 
posseses a nonnegative right equilibrium potential. Since f(U) is non singular we also 
have a left equilibrium potential, which we do not know if it is nonnegative. Then the 
first part is proven under the extra hypothesis that /(0) = 0. 

Assume now a = /(0) > 0, and consider g(x) = f(x) — a, which is a strictly increasing 
convex function. Obviously f(U) = g{U) + all', so 

1 1 

W> = TjTnF, ^(U) > 0, v f (u) = : f ff (lO, 

1 + a/i 9 ([/) 1 + ap g (u) 

where fi g (u)) = 1'a*s(J7)) > 0, and we have used the fact that p g (u)) — Vg(u))- Thus /(£/) 
has a nonnegative right equilibrium potential, and a left equilibrium potential. We need 
to prove that f(U) is nonsingular, and det(f(U)) > 0. This follows immediately from the 
equality 

f{U)=g{U){l + ap g{u) V). 

Indeed we have 

fill)- 1 = g(U)- 1 - p g{U ){v g {u))' and det{f(U)) = det(g{U))(l + ap g(p) ), 

1 + ap g (jj) 

from which the result is proven. 

In the bi-potential case use U' instead of U to obtain the existence of a nonnegative 
left equilibrium potential for f(U). 

□ 
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Proof. (Theorem 12.51) Using the same ideas as above we can assume that /(0) = 0. 
Also we have U~ l (l + tf(U)) = M t is a M-matrix, for all t > 0. Therefore I + tf(U) is 
nonsingular for all t, and we denote by fi t and v t the equilibrium potentials for I + tf(U). 

Assume first that f(U) is in class T which means that 



we obtain that \i t = M t 1 fi U > 0, because M t 1 is a nonnegative matrix. Thus, r(f(U)) = 
oo and f(U) is nonsingular. From Proposition II . II we get f(U) G biV. 

Reciprocally if f(U) G biV then r(f(U)) = oo and the result follows. □ 

Lemma 3.2 Assume that U G V. Then any principal square submatrix A of U is also 
in class V. The same is true if we replace V by biV. 

Proof. By induction and a suitable permutation is enough to prove the result for A the 
restriction of U to {1, . . . , n — 1} x {1, . . . , n — 1}, where n is the size of U. Assume that 



Since A~ x = A — we obtain that the off diagonal elements of A -1 are non-positive. 
It is quite easy to see that the result will follow as soon as A" 1 ! > 0. 
Since UePwe have that Al — ( > and 9 > g'l. Therefore, 



□ 

In what follows given a vector a we denote by D a the diagonal matrix, whose diagonal 
is a. 

Lemma 3.3 Assume U G biV and a > 1. If 



r(f(U)) 




We prove that for all t > 0, /i f , are nonnegative. Since 



M t fi t = U~ 1 1 = 




A" 1 ! = Al - -C^'l = Al - ^ C > Al - C > 0. 




then there exists a nonnegative vector 7] such that 



n 



b {a) . 
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Proof. We first perturb the matrix U to have a positive matrix. Consider e > and the 
positive matrix U e — U + ell'. It is direct to prove that 

= U- 1 - -J—^iyu)', 

where Jjlu = V \i v is the total mass of hjj. Then U e G biV and its equilibrium potentials 
are given by 

1 1 

l*U t = -i i — Vu, v Vt = , . — vu- 
1 + e/jLu 1 + ei/[/ 



We decompose the inverse of U e as 

ur 1 



0'e ^ 



and we notice that U e (e + #A = which implies that 

K = U e \ e , 

with X e = —j-Ce > 0. Also we mention here that X e is a sub-probability vector, that is 
l'A e < 1. This follows from the fact that U~ x is column diagonally dominant. 

Take now the matrix V e = D^ 1 U e . It is direct to check that V e G M.~ l and its 
equilibrium potentials are 

Hv e = K, v Vt = D be v Ue . 

Thus V t G biV and we can apply Theorem 12.11 to get that the matrix Ve posses a right 
equilibrium potential rj e > 0, that is, for all i 

£(K Ca) )«fo e )i = i, 

3 

which is equivalent to 
Hence 

U^T] e = b { e a) . 

Recall that the matrix U^ a > is nonsingular. Since obviously ut ] -> as e -> 0, we get 



and the result follows. □ 



Proof. (Theorem 12. 2L ) Consider first the case where U G biV. We already know 
that U {a) is nonsingular and that it has left and right nonnegative equilibrium potentials. 
Therefore, in order to prove that U {a) G biV is enough to prove that (?/(«)) -1 is a Z-matrix 
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that is for i ^ j we have ((U^) -1 )^ is non positive. An argument based on permutations 
shows that it is enough to prove the claim for % — l,j — n, where n is the size of U. 
Decompose and its inverse as follows 

We need to show that (3 > 0. We notice that 5 = > and that 

+ 5b^ = 0, 

which implies that 




Therefore, j — r\ > 0, where 77 is the vector given in Lemma 13.31 Thus [3 > and the 
result is proven for the case U e biV. 

Now consider U the inverse of the M-matrix M. Using Theorem 2.5.3 in jH], we get 
the existence of two positive diagonal matrices D, E such that DME is a strictly row 
and column diagonally dominant M-matrix. Thus V = E~ 1 UD~ 1 is in biV from where it 
follows that V {a) E biV. Hence, U (a) = E {a) V {a) D ia) is the inverse of an M-matrix. The 
rest of the result is proven in a similar way. □ 

Proof. (Theorem 12. 31 ) By hypothesis we have U = I — P, where P > and PI < 1. 

We notice that U is diagonally dominant on each row, that is for all i,j 

Also we notice that U = I + PU and therefore Uu > 1. 

According to Theorem 12.21 we know that H = (U^)" 1 is a row diagonally dominant 
M-matrix. The only thing left to prove is that the diagonal elements of H are dominated 
by one, that is Ha < 1 for all %. It is enough to prove this for % = n, where n is the size 
of U. 

Consider the following decompositions 




A direct computation gives that 

7 = P 5 + X 'P > pS. 
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Since by hypothesis 7 < 1, to conclude that 6 < 1 it is enough to prove that p > 1. We 
have 

p = (1 - Pnn )u: n - = tc - 5>n^ a n = - '• 

On the other hand we have U"~ l < U!^' 1 and J^PnjUjn = U nn — 1 from where we deduce 

j 

P > u^- 1 > 1. 

This finishes the case when PI < 1. The rest of the result is proven by using U' instead 
oiU. □ 



4 Proof of Theorem 12.4 



Let n be the dimension of U. Notice that U is a GUM if and only if n < 2 or every 
principal sub-matrix of size 3 is a GUM. 

Since by hypothesis the matrix I + tU is a bi-potential, it is diagonally dominant 

l + tU ti >tUij, 

and we deduce that Ua > Uij (by taking t — > oo). This proves the result when n < 2. So, 
in the sequel we assume n > 3. 

Consider A any principal sub- matrix of U, of size 3x3. Since I + tf(A) is a principal 
sub-matrix of I + tf(U), we deduce that I + tf(A) E biP (as long as I + tf(U) E UP). If 
the result holds for the 3x3 matrices, we deduce that A is GUM implying that U is also 
a GUM. 

Thus, in what follows we consider that U is a 3 x 3 matrix that verifies the hypothesis 
of the Theorem. After a suitable permutation we can further assume that 




where a = min{f4j : i ^ j} = min{[/} and (3 = min{Uji : Uij = a, i ^ j}. 

Since U is diagonally dominant we have min{a, d, e} > a. Take / increasing such that 
f(a) = and f(x) > for x > a. Then, 

'1 + f(a) f(h) f(b 2 ) 
l + f(U)= | f( Cl ) l + f(d) 

/(c 2 ) ftf) 1 + /(e) 



13 



is a biV-matrix whose inverse we denote by 

5 -pi -p 2 \ 
-di 7i -72 • 

-Q 2 _ 73 74 J 

In particular we obtain 

A + /(d) V 1 _ / 7l - 72 \ _ 1 /0A />A ' 
V /(/?) 1 + f(e)J ~ V-73 74^ 5 \e 2 J \ P2 ) ' 

and we deduce that 

= 72 = 01/*. (4-1) 



Case p 2 = 0. We deduce that /(&2) = 0, and then 

b 2 = a, and c 2 > f3, (4.2) 
where the last conclusion follows from the definition of (3. Therefore we have that 



U = ci d a , (4.3) 




and we should prove that U is GUM. 

Now consider another increasing function g such that g(/3) = and g{x) > for 
x > p. Then, 

'l + <7(a) 0(61) 
!+<?(£/)= I ^(ci) l + ^(d) 

o(c 2 ) 1 




Its inverse is of the form 




-Pi 

7i 
-7s 

As before we deduce that = 73 = 2 p 1 . 

— Subcase 9 2 = 0. We have g{c 2 ) = which implies c 2 = (3. Thus, in this 
situation we have that U is 
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By permuting rows and columns 1,2, if necessary, we can assume that b\ < c±. 
Consider the situation where c\ < j3, of course implicitly we should have a < (3. 
Under a suitable increasing transformation h we have 



l + h{U) 






1 + h(d) 
HP) 



and its inverse 



/ 



i 



l+h(a) 







1 



h(/3) 

\ (l+ft(a))(l+fc(e)) 



The sum of the third row is then 
1 



l+h{d) 

HP) 






-h(e) 



\ 





1 + h(e) 



1 - h{(3) 



(l+h(d))(l+h(e)) l+h(e)J 



+ 



1 + h(a) 1 + h(d) 



and this quantity can be made negative by choosing an appropriate function 
h. The idea is to make h(j3) — > oo and 



1. 



max{/i(a), h(d)} 
Therefore, a > (5 and U is a GUM. 

Subcase p\ = 0. We have g(bi) = and then b\ < (3. Take again an increasing 
function, denoted by £, such that 



i + e(u) 






1 + £{d) 
1 + £{e) 



and its inverse 




l 



{l+l(a))(l+i(d)) l+l(d) 

v e -kl o 

\ (l+e(a))(l+£(e)) U 





1 



l+£(e) 



The sum of the first column is 

1 



1 + £(a) 



1 



£( Cl ) 



[l+W) (l + i(e)) 



which can be made negative by repeating a similar argument as before, if both 
Ci > (3 and c 2 > (3. 
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So, if ci > j3 we conclude that c 2 < /3, but we know that c 2 > (3 (see ()4.2|) ) and 
we deduce that c 2 = f3. Thus, a < b\ < {3 < c\ and 



U 



which is a GUM. 

Therefore we can continue under the hypothesis c\ < (3 < c 2 . 
* Subsubcase b\ < (3. Again we must have a < (3. Under this conditions 




we have that c 2 > a. Using an increasing function k we get 
I + k(U) = 






1 + k(d) 
H(3) 




and its inverse is 
/ 



i 



l+k(a) 
fc(ci) 



1 



V 



(l+fc(o))(l+fc(d)) 

fc(c 2 )(l+fc(d))-fc(/3)fc(ei) 
'(l+fe(a))(l+fe(d))(l+fc(e)) 



l+k(d) 



The sum of the third row is 

1 (. k{c 2 ) 



(l + k(e)) 



1 



(l+fc(d))(l+fc(e)) l+fc(e) 



k{(3)k{ Cl ] 



K(3) 

1 + k(a) ^ (1 + k(a))(l + k{d)) 1 + k(d) 



If Ci < f3 we can assume that k(c\) = 0, and we can make this sum to be 
negative by choosing large k. Thus we must have c\ = (3, in which case 
the sum under study is proportional to 



1 - 



k(c 2 ) 



k(py 



HP) 



If c 2 = (3 then 



l + k(a) (1 + k(a))(l + k(d)) 1 + k(d) 



U 



(4.4) 




is a GUM. So we must analyze the case where c 2 > (3 in ()4.4j) . We will 
arrive to a contradiction by taking an asymptotic as before. Consider a 
fixed number A G (0, 1). Choose a family (/c r ) r gN such that as r — ► oo 



oo, 



fc r Q3) 

k r (c 2 ) 



A. 



fcr(c 2 ) 

/c r (a) 



k r (d) 
k r (a) 



1(3 



where = 1 if d > (3, and = A if d — (3. The asymptotic of (|4.4|) is then 

This quantity is strictly negative for the two possible values of 0, which is 
a contradiction, and therefore c 2 = (3. 

To finish with the Subcase p~\ = 0, which will in turn finish with Case 
P2 = 0, we consider 
* Subsubcase b\ = (3. We recall that we are under the restrictions c\ < 

(3 < c 2 and 

(a [3 a N 
C\ d a 
c 2 (3 e t 

Notice that if c 2 = f3 then U is GUM. So for the rest of this subcase we 
assume c 2 > (3. Also if c\ — a we can permute 1 and 2 to get 



mm' 



which is also in NBF, and U is a GUM. Thus we can assume that c\ > a, 
and again of course we have a < (3. 
Take an increasing function m such that 

/l + m(a) m(f3) 
I + m(U)=\ m{c x ) l + m(d) 

y m(c 2 ) m(P) 1 + m(e) 

We take the asymptotic under the following restrictions: 

MP) x c / n n , x m(e) m(c2) m(d) 

-— "> A £ (U, 1J, — r ► A, — — > 1, — r ► 1, -— — > <p, 

m{a) m{a) m(a) m{a) m{a) 

where = 1 if d > f3, and it is A if d — (3. The limiting matrix for 

V = 





whose determinant is A = — A 2 > 0. Therefore V must be in biV. On 
the other hand the inverse of V is given by 
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and the sum of the first column is 

A 2 -A 

which is a contradiction. 

This finishes with the subcase P2 = and we return to (|4.1|) to consider now the 
following case 

Case 6*i = 0. Under this condition we get c\ = a and 



a 


h 


b 2 


a 


d 


a 


("2 





e 



Consider the transpose of U and permute on it 2 and 3, to obtain the matrix 



U 



where now bi > (3. Clearly the matrix U verifies the hypothesis of the Theorem and 
has the shape of (|4.3|) . that is we are in the "case p2 = 0" which we already know 
implies that U is GUM. Therefore U itself is GUM. 



a 


C2 


a 


b 2 


e 


a 


bi 





d 



□ 



5 Filtered Matrices, sufficient conditions for classes 
biV and T 

A class of matrices of our interest is the class of filtered matrices, which turn to be a 
generalization of GUM matrices. They were introduced as operators in [2j to generalize 
the class of self adjoint operators whose spectral decomposition is written in terms of 
conditional expectations (see for instance 0, |lj and [TU]). 

The basic tool to construct these matrices are partitions of J n = {1, • • • ,n}. The 

components of a partition 1Z are called atoms. We denote by ~ the equivalence relation 
induced by 1Z. Then i, j are in the same atom of 1Z if and only if i ~ j. 

A partition 1Z is coarser or equal than Q if the atoms of Q are contained in the atoms 
of 1Z. We denote this (partial) order relation by 1Z =^ Q. For example in JZ4 we have 
1Z = {{1, 2}, {3, 4}} =4 Q = {{1}, {2}, {3, 4}}. The coarsest partition is the trivial one 
A/" = {Jn} and the finer one is the discrete partition T = {{!}, {2}, • • • , {n}}. 
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Definition 5.1 A filtration is an strictly increasing sequence of comparable partitions 
F = -< 72-1 -<•••-< TZk}- A filtration in wide sense is an increasing sequence of 
comparable partitions G = {TZq =4 TZi =4 ' ' ' TZk}- 

The difference between these two concepts is that in the latter repetition of partitions 
is allowed. 

Definition 5.2 A filtration ¥ = {TZq -<■■■-< TZk] is called dyadic if each non-trivial 
atom ofTZ s is divided into two atoms oflZ s+ \. 

The following example is the simplest dyadic filtration 

F = {M -< {{1}, {2, • • • , n}} -<■••-< {{1}, {2}, • • • , {*}, {i + 1, • • • , n}} ■< ■ ■ ■ ■< T}. 
Each partition 1Z induces an incidence matrix F =: F(TZ) given by 



F 



1 if % S j 
otherwise 



A vector v G M. n is said to be ^-measurable if v is constant on the atoms of TZ, that is 

. iz . 

I ~ J Vi = Vj . 

This can be expressed in terms of standard matrix operations as 

F(TZ)v = D wn v, 

where w-r, = F(1Z)1 is the vector of sizes of the atoms. Recall that D z is the diagonal 
matrix associated to the vector z. The set of 7?.-measurable vectors is a linear subspace 
of R n . 

Definition 5.3 A matrix U is said to be filtered if there exists a filtration in wide sense 
G = {<2o ^ Qi =4 ■ ■ ■ =4 Qi}, vectors ax, ■ ■ ■ , a/, b±, ■ ■ ■ , bi with the restriction that a s , b s 
are Q s+ i-measurable, such that 

i 

U = Y,D a F(Q s )D bs , (5.1) 

There is no loss of generality if we assume that Q = H and Qe = JF, that is F(Q ) = 11' 
and F(Qi) = I, the identity matrix. Let us see that (|5.1|) can be simply written in terms 
of a filtration. Indeed, notice that if a s and b s are Q s -measurable then 

D as F(Q s )D bs = D as D bs F(Q s ) = D asQbs F(Q s ), 
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where the vector a s 0b s is the Hadamard product of a s and b s which is also <2 s -measurable. 
Hence a sum of terms of the form 

D Ch F(Q s )D ba + D as+1 F(Q s+1 )D bs+1 + ■■■ + D as+r F(Q s+r )D bs+r , 

with TZ — Q s — ■ ■ ■ — Q s+r , can be reduced to the sum of two terms as 

D c F(K) + D as+ F(K)D bs+r , 

r-1 

where C = a s+ k b s+ k, which is 7£-measurable. In this way representation (j5.1J) can 

k=0 

be written as 

fc 

U = J2 D C F{K S ) + D m F(U s )D ns (5.2) 

where F = {TZo -< 7Zi -< ■ ■ ■ -< 7Zk} is a filtration, Af = 7Zo, T = TZk, C s is 1Z S - 
measurable, m s ,n s are 7^<j + i-measurable and = 0. We shall always consider this 
reduced representation of (|5.1jl . and we shall say that U is filtered with respect to the 
filtration F. 

If all m s ,n s are 7^,,-measurable then (|5.1|) reduces to the form 

k 

F(TZ S ), (5.3) 

s=0 

and U is a symmetric matrix. 

We are mainly interested in a decomposition like (|5.2j) with the vectors m s ,n s having 
the following special structure: 

m s = r s 0p s , n s = q s , (5.4) 

where T s is 7?. s -measurable and {p s ,q s } is a 7?. s+ i-measurable partition, that is, they are 
7£ s+ i-measurable {0, l}-valued vectors with disjoint support: p s q s = and p s + q s = 1. 
If this is the case we say that U is a Special Filtered Matrix (SFM) 

k 

U = Y,D Cs F{TZ S ) + D r D Ps F(K a ) D Qs . (5.5) 

Notice that = 0. 

It is not difficult to see that every CBF matrix is filtered. This is done by induction. 
Assume that 

TJ — ( ^ ^pln-p 

" \(31 n . p V p B 
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Define TZq — N and 1Z\ = {{1, • • ■ ,p}, {p + 1, ■ • • , n}}. Take Co = al n , T = (/? — 
a)l n ,p = (O p ,l n -p)',q = (l p ,O n -p)'. Then 

D Co F(K Q ) + D To D po F(n )D qo = ( * lpl \, n f pK v P ) ■ 

The key step is that A — a, B — a are also in CBF. We have that Co, T are 7^0-measurable 
and po, q is a 7^i-measurable partition. We also notice that if < a < [3 then C > 

o,r > o. 

k 

The induction also shows that U can be decomposed as in ()5.5|) U = ^2 Dc s F(lZ s ) + 

s=0 

Dr 3 D Ps F(7l s )D qg , where F = {TZq -<■■•-< IZk] is a dyadic filtration, C S ,T S are 1Z S - 
measurable and {p s , q s } is a T^s+i-measurable partition. 

We summarize now the representation form for the class of CBF, NBF and GUM 
matrices. 



Proposition 5.1 V is a permutation of a CBF if and only if there exists a dyadic filtra- 
tion F = {TZ -<•■■-< IZk}, o, sequence of vectors C , • • • , Ck, T , ■ ■ ■ ,T k verifying that 
for all i: C s , T s are lZ s -measurable and a sequence {p s , q s } of lZ s+ i-measurable partitions, 
such that 

k 

V = Y, D Cs F{K s ) + D rs D Ps F(TZ s )D qs , 

s=0 

that is V is a SFM. Also V is a permutation of an increasing CBF matrix if and only if 
there is a decomposition where T , C s , T s : s = 1 . . . , k, are nonnegative. On the top of 
this V is a nonnegative matrix if and only if Cq is nonnegative. 

Moreover, V is a GUM if and only if C s , T s : i — 0, • • • , k are non-negative, and for 
i = 0, . . . , k — 1 it holds 

r s <c s+1 + r s+1 . (5.6) 

Finally, V is an ultrametric matrix if and only if there is a decomposition with T s = for 
all s. □ 

Remark 5.1 We can assume without loss of generality that each p s ,q s is obtained as 
follows. The nontrivial atoms Ai, ■ ■ ■ A r oflZ s are divided into the new atoms 

"4-1,1) -4.1,2 " • • >Ar,l, Ar,2 

of 1Z s+ \. Consider also B\,---B r the set of trivial atoms in TZ S (that is the singleton 
atoms). Take q s be the indicator of Ai,i U • • ■ U Ar,i, Ps be the indicator of A± t 2 U • ■ ■ U 
A r ,2 U E\ U • ■ • U B r and T s = on B = B\ U • • • U B r , which is 1Z S -measurable. We point 
out that the partition 1Z s +i is obtained from 1Z S refined by p s . Also the important relation 
holds 

D Ps F(n s )p s = D Ps F{K s+1 )l (5.7) 
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Example 5.1 Consider the CBF matrix 



U 



fa CK2 Oi\ ot{\ 

@2 b Oil Oi\ 

Pi Pi c a 2 

\Pi Pi k dj 



U is a NBF if the following constraints are verified: a,\ < Pi, a.\ < mm{« 2 , 0:2}, Pi < 
min{P2, P2}, ct2 < P2, CX2 < P2 and finally the diagonal dominates over each row and 
column, that is P2 < min{a,b}, P2 < min{c,d}. 

U is filtered with respect to the dyadic filtration TZq = {1,2,3,4} -< TZi — {{1, 2}, {3, 4}} -< 
TZ 2 = {{1}, {2}, {3}, {4}} and can be written as 

U = D Co F(n ) + D To D po F(TZ )D qo + D^F^) + D^D^F^D^ + D C2 F(n 2 ), 



where 



( "A 

Oil 

\ ai / 
(a 2 — «i\ 



C'i 



a 2 

\&2 



Oi\ 
Oil) 



(Pi- a ± \ 
Pi - «i 

Pi ~ Oil 
\Pl -Oi X ) 

(P 2 - 0L 2 \ 
P 2 ~ 0L 2 
P 2 ~ 6i 2 

\(3 2 - a 2 J 



Po 



Pi 



/0\ 



1 

VV 

M 
1 


VV 



Qo 



Qi 



fl\ 
1 



W 

/A 


1 



and 



Co 



/ a 
b 

c 

\d 



a 2 \ 

0L2 

&2 

ol 2 ) 



The constrains are translated into: the positivity of these vectors and the ones induced by 
A5.6)) . We point out that we can also choose, for example, Ti = (0, P2 — 02, O,/^ — 0.2)' , but 
in this case T\ is not lZ\-measurable. As we will see in subsection \5.1\) this measurability 
condition will play an important role. 

Example 5.2 Consider the nonnegative CBF matrix 



U 



This matrix is a SFM and can be decomposed as in \5.d\) . Nevertheless no such decompo- 
sition can have all terms nonnegative. In particular no permutation of U is an increasing 
CBF. 
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Remark 5.2 Notice that the class of CBF is stable under Hadamard functions. Never- 
theless there are examples of filtered matrices for which f{U) is not filtered. Consider the 
matrix 

U = D a F 1 + D a F 1 D b + D (3 F 2 , 

where F\ = F(Af) = 11' and F 2 = I. We have a is a constant vector, and we confound 
it with the constant aGl The vectors a, b, (3 are all T ' -measurable. Then U is filtered 
and moreover 

U = a + ab' + Dp (5.8) 

Take a = (3 = and a = (2, 3, 5, 7)' and b = (11, 13, 17, 19)'. Then all the entries of 
U are different. As f runs all possible functions f(U) runs over all 4 x 4 matrices. This 
implies that some of them can not be written as in \5. because in this representation 
we have at most 13 free variables. Still is possible that each f(U) is decomposable as in 
\5. 1\) using maybe a different filtration. A more detailed analysis shows that this is not 
the case. For example if we choose the filtration Af -< {{1, 2}, {3, 4}} -< T then every 
matrix V filtered with respect to this filtration verifies that V13 = V23 = Vu = V24. 

Matrices of the type F(TZ) are related to conditional expectations (in probability the- 
ory). Indeed, let 1Z = {Ai,A2, • • • ,A r } and ng = #(Ae) be the size of each atom. It is 
direct that w = w-r = F(TZ)1 is a ^-measurable vector that verifies Wi = ne for i G At- 
Then 

E n = D~ l F(K) = F(K) D- 1 

is the matrix of conditional expectation with respect to the a-algebra generated by 1Z. 
This matrix E = E-^ satisfies: 

EE = E; E' = E; El = 1; 
Wv Kv is 71 — measurable; 
if v is 1Z — measurable, then Kv = v. 

Therefore, E is the orthogonal projection over the subspace of all 7£-measurable vectors. 
In the case of the trivial partition Af one gets E^/ - = -11' the mean operator. 

Remark 5.3 The L 2 space associated to {1, ■ • • , n} endowed with the counting measure, 
is identified with M. n with the standard euclidian scalar product. In this way each vector of 
W 1 can be seen as a function in L 2 , and E is an orthogonal projection. The product D V K 
(as matrices) is the product of the operators D v and E, where D v is the multiplication by 
the function v. Notice that KD V and E(-u) are quite different. The former is an operator 
(a matrix) and the latter is a function (vector). They are related by E(-u) = E,D V (1), 
where 1 is the constant function. 
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Let 7Z, Q be two partitions, then 1Z ^ Q is equivalent to E-^Eg = EgE^ = E-#. This 
commutation relation can be written as a commutation relation for F(TZ) and F(Q). In 
fact, 

F(K)F(Q) = E n D Wn E Q D WQ = E n E Q D WK D WQ = 

= E n D Wn D WQ =F(n)D WQ , 
F{Q)F{K) =(F(n)F(Q))' = D WQ F(Tl). 



5.1 An algorithm for filtered matrices: conditions to be in biV 

In this section we explain a backward algorithm to determine when a filtered matrix is in 
class biV. Assume that U has a representation as in (|5.1j) 

£ 

s=0 

where we assume further that a s , b s are all non-negative. In particular U is a non-negative 
matrix. 

We introduce the conditional expectations E s = E Qs = D^Lj F(Q S ) and the nor- 
malized factors: a s — a s © F(Q S )1, b s = b s . Then U can be written as 

£ i 
U = Y,D aa E s D ba = ^a s EA, 

s=0 s=0 

where we have identified vectors (functions) and the operator of multiplication they in- 
duce. We shall use this notation throughout this section. Finally, we remind that E^ = I. 

We can now use the algorithm developed in [7] to study the inverse of I + U. In 
what follows, we take the convention 0-oo = 0/0 = 0. This algorithm is defined by the 
backward recursion starting with the values A^ = fig — Kg — 1, og — (1 + agbg)~ x and for 
s = 1- I,- - ,0: 

A s = A s+ i[l — cr s+ ia s+ iE s+ i(hi s+1 b s+ i)]] 

= /^s+ifl — (?s+ib s +iE s+ i(K s+ ia s+ i)]; (^q) 
k s = E s+ i(X s ) = E s+ i(/i s ); 

a s = (1 + E s (n s a s b s )y 1 , 

from where we obtain the recursion 

^ , , EjKM H )EjK H bc) 

= EJk s ) - sK s % sK 7 S N . 5.10 
1 ; l+E s (K s a s b s ) 1 ; 

The algorithm continues until some A or fi are negative otherwise we arrive to s = 0. If 
this is the case then I + U is nonsingular and its inverse is of the form I — N where 

£ 

N = 22 o" s A s a s E s 6 s /x s . 
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We also have that 

A_i = (I — N)l, and = (I — N)'l, 

where A_i,/i_i are obtained from the first two formulae in (|5.9|) for s = —1. Therefore, if 
they are also non-negative the matrix I + U is a biV-m&trix. 

In this way we have that a sufficient condition for 1 + U to be a foP-matrix, is that the 
algorithm works for s = £,■■■, and all the A, /i are nonnegative, including A_i, \x-\. In 
this situation we have that A (and \i) is a decreasing non-negative sequence of vectors. 

Sufficient treatable conditions involve the recurrence ()5.10|) . Starting from Kg = 1 we 
assume this recurrence has a solution such that k s G [0, 1] for all s = £, ■ • ■ , — 1. We shall 
study closely this recursion for the class of SFM, and we shall obtain sufficient conditions 
to have I + U in biV. 

Before studying this problem, we discuss further the algorithm. We have the following 
relations: 

( t \ x i 

( I + Yu a Ah J = I - Z crkhaAhHk = I - N s , 

\ k=s / k=s 

A s _x = {I-N s )l, // s _! = (I-N s )'l 

That is, the algorithm imposes that all the matrices: 

i i 
I + aiEebe, • ■ ■ , I + ^ a k E k b k , ■ ■ ■ , I + ^ a k E k b k —I + U, 

k=s k=0 

are in class biV. 

We now assume that U is a SFM with a decomposition like 

k 

U = Y,D Cs F{K S ) + D Vs D Ps F(TZ S ) D qs , 

s=0 

where F = TZq -<■■■-< lZ k is a filtration, C S ,T S are nonnegative 7£ s -measurable and 
{p s ,Qs} is a 7£ s+ i-measurable partition. Again we put E s = D~^, n ^ x F(1Z S ) and the 
normalized factors: 

c s = C s F(K.)1, is = T s F(TZ S )1, 

which are 7£ s -measurable. Since diagonal matrices commute we get that U has a repre- 
sentation of the form 

k 

U = c s E s + 7 s p s E s g s , 

s=0 

with 7^ = 0. In the previous algorithm we can make two steps at each time and consider 
k s in place of k 2s , X s instead of A2 S +i, l s instead of A2 S . We also introduce d s = 1/k s 
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to simplify certain formulae (this vector can take the value oo). We get, starting from 
Kk — h — 1) Cfc = (1 + Cfc)~\ that for s — k — 1, • ■ • ,0 

A s = O's+lh+l] 

l s = A s [i - jsPsEs^s/ (c s+ i + 4+1))]; 

k s = E S (Z S ); 

cr s = 1/(1 + k s c s ) = d s /(c s + d s ). 

Similar recursions hold for /i, m, which are the analogous of A, I. Relation ()5.10|) takes the 
form 

t = Es ( — TT-) ~ 7 * E * ( — h~) E * ( — h~ ) ( 5 - n ) 

d s V c -s+i + d s+ i J V c s+ i + d s+ i J \ c s+ i + d s+ i J 

The inverse of I + U is I — N where 

k k-l k 

N = ^2 c s cr s l s E s m s + ^ j s X s p s E s q s fi s = ^ c s a s l s E s m s + j s X s p s E s q s fi s . (5.12) 

s=0 s=0 s=0 

Again A_i = (I — N)l = s l , and similarly = s m . 
Let us introduce the following function 

p s = E s (p s )p s + E s (q s )q s . 

Theorem 5.1 Assume that the backward recursion H5.ll)) has a non-negative solution 
starting with df. — 1. Assume moreover that this solution verifies for s = k — 1, • • • ,0 

Psls < c s +i + d s+ i. (5.13) 

Then \ s , l s , /i s , m s , s s : % — k, ■ ■ ■ , 0, as well as A_i, are well defined and nonnegative. 
Therefore, I + U G biV and its inverse is I — N where N is given by \5.1fy) . 

The proof of this result is based on the following lemma. 

Lemma 5.1 Assume x,y are nonnegative vectors and E is a conditional expectation. If 
xE(y) < 1 then E(xy) < 1. 

Proof. We first assume that y is strictly positive. Since x < 1/E(y) and E is an increas- 
ing operator, we have 

For the general case consider (?/ + el)/(l + e|x| 00 ) instead of y and pass to the limit e — > 0. 
□ 

Proof. (Theorem 15. 1|) We notice that condition ()5.13|) implies that 

-E s (q s ) < 1. 



c s +i + d s+ \ 
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Since 7 S is E s -measurable and q s = q 2 s we obtain 

^(—h-)=M^r-)- 

\c s+ i + d s+ ij \c s+ i + d s+ ij 

This last quantity is bounded by one by Lemma f5. II Similarly we have 

( — x^— J ^ X > 

\ c s+l ~r a s+l J 

which implies that the algorithm is not stopped, all the coefficients are non-negative 
including A_i, □ 



Corollary 5.1 Assume that for s — k — 1, • • • , we have 

Psls < c s+ i + 7s+i- ( 5 - 14 ) 

Then the recursion 115. 11]) has a nonnegative solution that verifies h5.1ty . In particular, 
I + tU is in class biV for all t > 0, and U is in biV if it is nonsingular. 

Proof. Let us consider first the case t = 1. We prove by induction that 7 S < d s . For 
i = k we have = 7^ < dk — 1. We point out that if we multiply in (|5.11|) by 7 S we get 



d s \c s +i + d s+ ij \ c s+i + d s+ i J \c s +i + d s+ i 

which is of the form x + y — xy, where x = E s ( — Y The inequality f)5.14|) . the 

induction hypothesis 7 s+ i < d s+ ± and Lemma 15.11 imply 0<x<l,0<y<l. In 
particular 

< ?r < 1, 
d s 

and the induction is completed. Theorem 15.11 shows that I + U is in class biV. We notice 
that tU also verifies condition (J5.14J1 because this condition is homogeneous, and the result 
follows. □ 

Remark 5.4 We notice that condition \5.14 ) can be expressed in terms of the original 
coefficients C,F in the dyadic case. In fact (see \5. r /\ ) ) 

p s ~E s (p s ) = DpiD'ytji F(n s )p s = D Ps D F l {ns)1 F(n s+1 )l, 

which implies that 

Ps = (l/F(K s )l)QF(K s+1 )l. 

Then inequality \5.14 ) is 

< c. s+ i + r s+ i, 

which is the condition for having a GUM (see \5.b)) ) . We mention here that condition 
\5.1J$ is more general than having a GUM, as the following example shows. 
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Remark 5.5 Consider the matrix Up 

(\ 0\ 
10 

p 13 1 
V? /? 1/ 

where TZ = Af,T = (3(1, 1, 1, 1)' < Ci = (1, 1, 1, 1)'. We compute c = 0, 70 = 4/?, a = 
Ci, 71 = and also po = 1/2. 

direct to check that Up 1 = U-p. Then for all (3 > i/ie matrix Up G M~ x . Also 
Up G foP and onfo/ i/0 < /5 < 1/2. WTien /5 > £/&e condition 115. b)) . that is 

is equivalent to (3 < 1. Then, this condition does not ensure that U G ozT 7 /i/ws happens 
because the filtration is not dyadic). Nevertheless, the analogue condition in terms of the 
normalized factors \5.14 ) 

Polo < ci +71, 
which is equivalent to (3 < 1/2, is £/ie ng/ii one. 

Corollary 5.2 Assume that 

k 

r=s+l 

hold for s = k — 1, ■ ■ ■ ,0. Then the recursion H5.ll)) has a nonnegative solution that 
verifies \5.1ty) . In particular, I + tU is in class biV for all t > 0, and U is in biV if it is 
nonsingular. 

Proof. Consider the set of inequalities 

Psls V £ s < c s+1 + £ a+1 , 

for i = k—1, ■ ■ ■ ,0. A non-negative solution is given by £ s = sup{0, ioPo~ Sr=i Cr > ' ' " ' IkPk 
Ylr=k+i c ri ' ' ' 1 1s-iPs-i — c s} ■ The hypothesis is that ^ = 0. Moreover we have that £ s 
is 7^-measurable. 

We show, using a backward recursion that £ s < d s . Indeed, l/£ s = E s (l/£ s ) > (c s+1 + 
^ s+ i) _1 by construction while l/d s < E s ((c s +i + d s+ i) _1 ). Then the inequality p s j s < 
Cs+i + Cs+i implies p s i s < c s+ i + d s+ i, from where the result holds (see Theorem 15. lj) . □ 



5.2 Conditions for class T and proof of Theorem 12.61 

Theorem 5.2 Assume that U has a decomposition 

1 

U = ^2a s E s b s , 

s=0 
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where a s , b s are nonnegative E s+ i-measurable. Then U belongs to the class T and more- 
over 

t(U) = inf{t > : (I + tUyH ^ or l'(I + til)' 1 ? 0} 
In particular if t(U) < oo then I + r(U) U G biV. 

Remark 5.6 Since the set of nonsingular matrices is open, then in the previous result 
when t(U) < oo, we have for t > t(U) sufficiently close to t{U), that the matrix 1 + tU 
is nonsingular. 

Theorem 15.21 states that every filtered matrix, with a nonnegative decomposition, is in 
class T which proves Theorem 12.61 

Proof. (Theorem 15. 21 ) A warning about the use of vectors and functions. Here 
we consider vectors or functions on {1, . . . ,n} indistinctively. Thus for two vectors a, b 
the product ab makes sense as the product of two functions, which corresponds to the 
Hadamard product of the vectors. Also an expression as (1 + ab)~ 1 is the vector whose 
components are the reciprocals of the components of 1 + ab. We also recall that (a)j is 
the 2-th component of a. 

First, for p = 0, . . . , i consider the matrices 

i 

U(p) = J2 a s E s b s . 

s=p 

We notice that U(0) = U . We shall prove that t p = r(U(p)) is increasing in p and t# = oo. 

We rewrite the algorithm for I + tU . This takes the form Xe(t) = fie(t) = Ke(t) = 
1, cre(t) — (1 + 1 a(b()~ l and for p — I — 1, • • • ,0: 

\{t) = 
H P (t) = 

K p (t) = 

a p (t) = 

Also A_i(t), /i-i(t) are defined similarly. If X s (t), fi s {t), cr s (t) : s — £,..., p are well 
defined then 

(I + tUip))- 1 = I-N(p,t), 

where 

i 

N(p,t) =J2°s(t)\ s (t)ta s E s b s Vs(t). (5.17) 

s=p 

If X s (t), fi s (t), a s (t) : s — £,... ,p are nonnegative then N(p,t) > 0, and (I + tU{p)) G 
M.~ l . Moreover, A p _i(t) and yU p _i(t) are the right and left equilibrium potentials of 
(I + tU(p)) 

(I + tU(p))X p ^(t) = 1, and + tU(p)) = 1'. 



X p+1 (t)[l - ap +1 (t)tap +1 E p+1 {K p+1 (t)bp + x)\; 
fj, p+1 (t)[l - o-p +1 (t)tbp +1 E p+1 (Kp +1 (t)a p+1 )}] 
E p+1 (X p (t))=E p+1 (ii P (t)); 
(l + E p (KJt)ta p b p ))-\ 



(5.16) 
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So, if they are nonnegative, we have I + tU(p) G biP. In particular we have that 

(I + taiEibe)- 1 = (l + tU(t))~ l = I - t(l + 1 a*&<) _1 a* E t b e . 

Since = I we obtain that Xg^i = \it-\ = (l + ta^) -1 . This means that I + tU(£) G biV 
for all t > 0. Therefore ti = oo and the result is true for U(£). This implies in particular 
that r^_i < T£. Assume the following inductive hypothesis 

• Tp+i < ■ ■ - < n ; 

and for q = p + 1, . . . ,£ 

• r g = inf{t > 0: Xq-i(t) ^ or /x,_i(t) ^ 0} = inf{t > 0: A g _i(t) ^ or ^ 0}; 

• A s (t), /x s (t), for s = £, . . . , q — 1, are strictly positive for t G [0, r q ) ; 

• If T q < oo we have I + r 9 {7 (g) G fezP. 

The case r p+1 = oo is clear. Indeed, fix t > 0. From Lemma l3~Tl 1 + tU(p+ 1) G 6iP 

and its equilibrium potential are strictly positive, that is X p (t) > 0,/i p (t) > 0. Thus, 

I + tU(p) is nonsingular, its inverse is I — N(p,t), where N(p,t) > is given by (|5.17|) . 
Hence, I + tU(p) G M^ 1 . We conclude that 

t p = inf{t > : I + tU(p) £ biV} = inf{t > : A p _i(t) ^ or //p_i(t) ^ 0}. 

So, if t p = oo we have, from Lemma f3. II that 

X p -i(t) > 0, fip-^t) > 0, 

and the induction step holds in this case. 

Now if t p < oo, by continuity we have I + T p U{p) G biV, and we shall prove later on 
that Ap_i(i), /x p _i(£) are strictly positive in [0,r p ). 

We analyze now the case r p+ i < oo. We first notice that in the algorithm the only 
possible problem is with the definition of a p (t). Since cr p (r p+ i) > the algorithm is well 
defined, by continuity, for steps £, . . . ,p on an interval [0, r p+ i + e], for small enough e. 
This proves that the matrix I + tU (p) is nonsingular in that interval, and that A p _i, yU p _i 
exist in the same interval. 

Now, for a sequence t n [ r p+1 either X p (t n ) or /i p (t n ) has a negative component. Since 
there are a finite number of components we can assume with no lost of generality that for a 
fixed component i we have (X p (t n ))i < 0. Then, by continuity we get that (A p (r p+1 ))j = 0, 
which implies (by the algorithm) that (A p _i(r p+1 ))j = 0. 

Assume now that for some t > r p+1 the matrix I+tU(p) G biV. Again by Lemma l3~T1 we 
will have that I + T p+1 £/(p) G biV but its equilibrium potential will satisfy A p _i(t p+1 ) > 0, 
which is a contradiction. Therefore we conclude that r p < r p+1 . 
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The conclusion of this discussion is that the matrix I + tU(p), for t G [0,r p +i], is 
nonsingular and its inverse is I — N(p, t), with N(p, t) > 0. That is I + tU(p) G Ai^ 1 and 
therefore 

r p = inf{t > : I + tU{p) £ biV} = M{t > : A p _i(t) ^ or ,u p _i(f) ^ 0}, 

and by continuity I + r p U (p) G biV. 

To finish the proof we need to show that r p coincides with 

S = inf{t > : Vi(*) > or ^ °}- 

It is clear that S < r p . If S < r p then, due to Lemma l3~Tt we have that both X p ^i(S) > 
and fip-i(S) > 0, which is a contradiction and then S = r p . This shows that A p _i(£), /x p _i(t) 
are strictly positive for t G [0, r p ), and the induction is proven. □ 



Remark 5.7 It zs possible to prove that k p (t p ) > w/ien r p < oo, but this is not central 
to our discussion. 
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